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Abstract 

It is due to the modularity of the analysis that results for cascaded systems have proved their util¬ 
ity in numerous control applications as well as in the development of general control techniques based 
on “adding integrators”. Nevertheless, the standing assumptions in most of the present literature 
on cascaded systems is that, when decoupled, the subsystems constituting the cascade are uniformly 
globally asymptotically stable (UGAS). Hence existing results fail in the more general case when the 
subsystems are uniformly semiglobally practically asymptotically stable (USPAS). This situation is 
often encountered in control practice, e.g., in control of physical systems with external perturbations, 
measurement noise, unmodelled dynamics, etc. This paper generalizes previous results for cascades 
by establishing that, under a uniform boundedness condition, the cascade of two USPAS systems 
remains USPAS. An analogous result can be derived for USAS systems in cascade. Furthermore, 
we show the utility of our results in the PID control of mechanical systems affected by unknown 
non-dissipative forces and considering the dynamics of the DC motors. 


1 Introduction 

Cascaded dynamical systems appear in many control applications whether naturally or intentionally 
due to control design. Cascades-based control consists in designing the control law so that the closed 
loop system has a cascaded structure. Such strategy has often the advantage of reducing the complexity 
of the controller and the difficulty of the stability analysis (see e.g. PEiil) as opposed to more general 
Lyapunov-based control methods. From a theoretical viewpoint the problem of stability analysis of 
cascaded systems has attracted the interest of the community starting with the seminal paper 23j. See 
also m and references therein. In general terms, a fundamental result that one may retain from the 
literature is that cascades of uniformly globally asymptotically stable systems (UGAS) remain UGAS 
if and only if the solutions are uniformly globally bounded (UGB). See [HI El] for the proof of this 
statement in the case of autonomous systems and m for the case of time-varying systems. In a 
similar spirit, see m for a local result: a proof that a LAS system perturbed by a converging input 
which is such that its solutions remain in the domain of attraction of the nominal system, remains 
LAS. Other works on stability of cascades deal, directly or indirectly, with the fundamental problem 
of establishing conditions for (uniform global) boundedness of the trajectories. A sufficient condition is 
the well-understood notion of input-to-state stability (ISS). 

A considerable drawback of most results on stability of cascades available in the literature is that 
they rely on the assumption of global stability properties for the separate subsystems. Nevertheless, 
in practice, it is often the case that only local (with a specified estimate of the domain of attraction) 
or semiglobal properties can be concluded. Semiglobal asymptotic stability pertains to the case when 
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one can prove that, by tuning certain parameter of the control system, the estimate of the domain of 
attraction can be arbitrarily enlarged. Such parameter is often, but not always, a control gain. We 
are not aware of results on semiglobal asymptotic stability nor ISS for cascades of “semiglobal ISS” 
systems. Another case which is not covered by most results in the literature of cascaded systems is 
that of stability with respect to balls. A notable exception, in a way, is [3] which introduced the notion 
of input to state practical stability. However, this notion is different from the one used here in the 
sense that the neighborhood of the origin which is stable is not required to be arbitrarily reducible by 
a convenient choice of a parameter. More generally, we are not aware of results on semiglobal practical 
asymptotic stability of cascades for continuous time-varying systems. For (parameterized) discrete-time 
systems the only results that we know of are those presented in m- Needless to say that the nature of 
discrete-time non-autonomous systems parameterized in the sampling time is fundamentally different. 

In this paper we address the stability analysis problem for cascades of time-varying systems that are 
uniformly semiglobally practically asymptotically stable (USPAS). We establish that, under a uniform 
boundedness condition on its solutions and provided a Lyapunov function for the “perturbed” subsys¬ 
tem, the cascade of two USPAS systems remains USPAS. In the same way, the cascade of two USAS 
systems is USAS. Our main result extends in this direction, HSJ Lemma 2] and the main results of 
HH1I2I! which are, in a way, at the basis of many theorems on UGAS of cascades. 

The rest of the paper is organized as follows. In next section we present some definitions of stability 
and an auxiliary proposition on semiglobal practical asymptotic stability. Our main result is presented 
in Section [SI In Section 0 we illustrate the utility of our findings with an example dealing with the PID 
control of a mechanical system taking into account the dynamics of the DC motors. The proofs of all 
the results are given in Section [5| and we conclude with some remarks in Sectional 


2 Definitions and preliminary results 


Notation. A continuous function a : M>o —► R>o is of class /C (a E JC), if it is strictly increasing and 
a(0) = 0; a E /Coo if) in addition, a(s) —> oo as s —> oo. A continuous function a : M>o —> M>o is of 
class C (cr E C) if it is non-increasing and tends to zero as its argument tends to infinity. A function 
(3 : M>o x M>o —► M>o is said to be a class ICC function if, /?(-,/) E /C for any t > 0, and /3(s, ■) E C 
for any s > 0. We denote by x(-,to,xo) the solutions of the differential equation x = f(t,x ) with initial 
conditions (t o,xo). We use |-| for the Euclidean norm of vectors and the induced L 2 norm of matrices. 
We denote by t3§ the closed ball in M n of radius 5, i.e. Bs := {x E M. n | |x| < <5}. We use the notation 
7i(5, A) := {x E | <5 < |x| < A}. We define |x| 5 := \x — z\. When the context is sufficiently 

explicit, we may omit to write the arguments of a function. 


2.1 Asymptotic stability of balls 

For nonlinear time-varying systems of the form 

x = f(t,x), (1) 

where x E M n , t E M>o and / : M>o x —> M n is piecewise continuous in t and locally Lipschitz in x, 

we introduce the following. 

Definition 1 (US of a ball) Let 5 and A be nonnegative numbers such that A > <5. The ball Bs is said 
to be Uniformly Stable on £>a for the system © if there exists a class /Coo function 77 such that the 
solutions of © from any initial state xq E £>a and initial time to E M>o satisfy 

\x(t,t 0 ,x 0 )\ s <v(\x 0 \), Vf>f 0 - 

Definition 2 (UA of a ball) Let 6 and A be nonnegative numbers such that A >5. The ball Bs is said 
to be Uniformly Attractive on B/\ for the system © if there exists a class C function a such that the 
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solutions of © from any initial state xq G Ba and initial time to G M>o satisfy 

\x(t,t 0 ,x 0 )\ 5 < a(t - t 0 ), Vt>t 0 . 

Definition 3 (UAS of a ball) Let 5 and A be nonnegative numbers such that A > <5. The ball 13 $ is said 
to be Uniformly Asymptotically Stable on Ba for the system © if it is both US and UA on Ba- 

Remark 1 The property of UAS of a ball defined above is less restrictive than the time-varying adap¬ 
tation of “asymptotic stability with respect to a set” given in [Hj for the case when the set is a ball. 
Indeed, in the latter reference, it is imposed that the ball B$ be positively invariant. Notice also that, 
modulo that Ba and B$ are closed, Definition 01 is equivalent to A —> 5 stability as defined in ©j]; more 
precisely, UAS of B$ on Ba implies A —> <5 stability and A —► 5 stability implies UAS of B$ on Ba>, for 
all A' G (6, A). 

We also emphasize that the UAS of a ball can be characterized by a ICC bound as in the case of UAS 
of the origin. 

Proposition 1 (ICC estimate for UAS of a ball) The ball B$ is UAS on Ba for the system HI) if and only 
if there exists a class ICC function (3 such that the solutions of © from any initial state xo G Ba and 
initial time to £ R>o satisfy 

\x(t,t 0 ,x 0 )\ s < (3(\x 0 \,t-t 0 ) , Vt > t 0 • 

Our main result for cascades is formulated based on the following property of boundedness of solutions. 

Definition 4 (UB) The solutions of © are said to be Uniformly Bounded on the compact set A C M n 
if there exist a class 1C function 7 and a nonnegative constant /r such that, for any initial time to G M>o 
and any initial state xo G A, it holds that 

[x(t,t 0 ,x 0 )| < 7 (|x 0 |) +/x, Vt > t 0 • 

2.2 Semiglobal practical asymptotic stability 

Our main result addresses the problem of uniform semiglobal practical asymptotic stability (USPAS) 
for parameterized nonlinear time-varying systems of the form 

x = f(t,x,9), ( 2 ) 

where x G M n , t G M>o, 9 G is a constant parameter and / : M>o x R n x W n —> M n is locally 
Lipschitz in x and piecewise continuous in t. 

Definition 5 (USPAS) Let 0 C M m be a set of parameters. The system © is said to be Uniformly 
Semiglobally Practically Asymptotically Stable on 0 if, given any A > <5 > 0, there exists 9* G 0 such 
that B$ is UAS on Ba for the system x = f(t, x, 9 *). 

In view of Proposition^ the above definition is equivalent to the following statement: for any given 
pair A > 5 > 0, there exists a parameter 9* G 0 and a ICC function f3 such that, for all xo G Ba and all 
to G M>o, | x(t, to, xo, 9*)\$ < /3(|xo| 3 1 —to) for all t > to- We stress that, although this does not explicitly 
appear in the notation, the function (3 is not required to be independent of 5 and A, as opposed to 
some other definitions of semiglobal and/or practical stability existing in the literature, such as [ 251 , Ej. 
This non-uniformity in 5 and A makes our definition of USPAS a more general property since, as it will 
appear more clear in the sequel, it allows to make use of Lyapunov functions with bounds that may 
depend on the tuning parameter 9. 

Systems of the form © result, for instance, from control systems in closed loop ; in this case, we 
can think of 9 as control gains or other design parameters. Then we say that © is USPAS if the 


estimate of the domain of attraction B/\ and the ball B§ which is UAS can be arbitrarily enlarged 
and diminished, respectively, by a convenient choice of the design parameters. Such situation is fairly 
common in control practice; for instance, in output feedback tracking control of mechanical systems 
(c/. m)- For discrete-time systems one also finds it useful to define USPAS with respect to a design 
parameter: in this case, the sampling time. See mm for definitions and a solid framework on USPAS 
for discrete-time systems. 

Notice also that when, by an abuse of notation, 5 = 0 we recover from Definition 0 the notion of 
uniform semiglobal asymptotic stability (USAS). If, in addition, A = oo, we recover the definition of 
uniform global asymptotic stability (UGAS). 

The following result gives a sufficient condition, in terms of a Lyapunov function, for the dynamical 
parameterized system © to be uniformly semiglobally practically asymptotically stable on a given set 
of parameters. See Section for the proof. 

Proposition 2 (Lyapunov sufficient condition for USPAS) Suppose that, given any A > 5 > 0, there 
exist a parameter 6 G 0, a continuously differentiable Lyapunov function 1 V : M>o x HU 1 —> M >0 , class 
/Cqo functions a 5 A , as, a, ay, a> and a continuous positive nondecreasing function c$ t a such that, for all 
x G 7i(6, A) and all t G M>o, 

< V(t,x) < 3i iA (|x|) (3) 

~Bt ^'^ + X ^ (yt ' x ’^ ~ “ a<5 ’ A (l x D ( 4 ) 

Assume further that, for any A > 0, 

limo^ oa g , A (S) = 0 (5) 

5 —>0 ’ 

and that, for any 5 > 0, 

lim a^o a s a (A) = oo . (6) 

A—>oo ’ 

Then, the system 0 is USPAS on the parameter set 0. 

It is worth mentioning that the condition a 5A (|®|) < V(t,x ) < ayAd^l); which implies ©, often 
holds in the analysis of control systems. In particular, it holds for systems with additive bounded 
disturbances when USPAS may be inferred using a Lyapunov function for UGAS of the corresponding 
unperturbed system. See Section 0 for an example. Condition © also appears naturally in the context 
of stability of perturbed systems. The last two conditions, © and ©, need to be imposed due to the 
fact that the Lyapunov function V is not required to be the same for all 5 and all A. As we show in 
the proof, conditions © and © “compensate” this lack of uniformity to ensure that the estimate of 
the domain of attraction B/\ and the set B$ which is UAS can be arbitrarily enlarged and diminished 
respectively. 

Remark 2 By noticing that the UAS of Bs on £>a implies the UAS of Bg’ on £>a' for any 5' and A' 
satisfying 5 < 5' < A' < A, the conclusion of Proposition El remains valid if © and 0 hold for all 5 
small enough and all A large enough. This relaxed assumption may be useful in practice. 


3 Stability of cascades 

We now consider cascaded systems of the form 

xi = fi(t,xi,9i) +g(t,x,0)x 2 (7a) 

X 2 = f2(t,x 2 ,0 2 ) (7b) 

x It should be clear that V may depend on S and A as well. We shorten the notation Vs ,a to just V for clarity. 
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where x := (x\ , x 2 ) G M ni x M™ 2 , 9 := (6j , 9 2 ) G M ni x M m2 , i G M>o, fi, /2 and g are locally Lipschitz 
in the state and piecewise continuous in the time. In order to simplify the statement of our main result, 
we first introduce the following notation. 

Definition 6 (V-set) For any A > 5 > 0, the V-set of 0 is defined as 

V f (S. , A) := {9 G M m : B s is UAS on B A for 0} . 

Theorem 1 Under Assumptions 00 below, the cascaded system 0 is USPAS on ©i x © 2 . 

Assumption 1 The system (I7bl) is USPAS on © 2 - 


Assumption 2 (Lyapunov USPAS of the xi -subsystem) Given any Ai > <5i > 0, there exist a parameter 
0^(<5i, Ai) G ©i, a continuously differentiable Lyapunov function Vi, class /Coo functions a$ 1 A , dj^Ai, 
Ai and a continuous positive nondecreasing function c^Ai such that, for all xi G TL{8\, Ai) and all 
t G M>o, bounds ©, 0 and 


8V '„ , 
9 ^ U ’ Xl) 


< CSi.Aid®!!) 


( 8 ) 


hold and, for any positive Ai, 

hin a^ Ai oa JliAl (^i) = 0. 


(9) 


Assumption 3 (Boundedness of the interconnection term) The function g is uniformly bounded both in 
time and in parameter, i.e. there exists a nondecreasing function G such that, for all x G M ni x M n2 , 
all 9 G ©1 x ©2 and all t G M>o, 

\g(t, x, 9) | < G{\x\). 


Assumption 4 (Boundedness of solutions) There exists a positive constant Ao such that, for any given 
positive numbers <5i, Ai, 5 2 , A 2 , satisfying Ai > max{<5i;Ao} and A 2 > <52, and for the parameter 
#*(<5i, Ai) as defined in Assumption^ there exists a parameter 9\ G Vf 2 (52, A 2 ) D @2 (cf. Definition^ 
and a continuous function 2 7 A 1 ,A 2 : ^>0 x ^>0 —> K>o such that 

lim 7 Ai ,a 2 (Ai, A 2 ) = +00, (10) 

Ai ,A 2 —>00 

and the trajectories of 0 with 9 = 9 * satisfy 

|®o| < 7 Ai,a 2 (Ai, A 2 ) =>- \x(t,t o ,x o ,0*)\ < Ai, Vt>t 0 . 

The proof of Theorem consists in constructing the balls Bs and B A and a KLC estimate for the 
solutions of the cascaded system, based on the respective balls for the x\ (i.e. (17al) with x 2 = 0) and 
the x 2 subsystems. For clarity of exposition we present this in Section f5. II 

In view of Proposition [21 Assumption [21 corresponds to the Lyapunov sufficient condition for USPAS 
of the zero-input xi-subsystem, with the additional condition of a bound on the gradient of Vi; we 
stress that the requirement corresponding to is no longer needed under Assumption 01 We state the 
main result under the more restrictive assumption than simply “USPAS” since our proof relies on the 
explicit knowledge of the Lyapunov function Vj. Besides J2J , we are not aware of a converse theorem for 
USPAS, as defined here, that gives all the required properties. In this respect, a converse theorem for 
USPAS of time-varying systems follows from m Corollary 1] for the case of locally Lipschitz /j, but 
it only establishes the two first inequalities of Assumption [21 and does not establish that the functions 
used to bound Vj can be chosen such that @ holds. Yet, both 0 and © are little restrictive, and are 
satisfied in many concrete applications, as for instance the case study of Section 0 

2 It should be clear that 7 may depend on <$1 and 62 as well, but we do not write this dependency explicitly in order to 
lighten the notation. 





Remark 3 In view of Remark [21 it is in fact sufficient that the requirements of Assumption [21 hold for 
all small 5 and all large A. Also, it is worth pointing out that Assumption ^ may be relaxed to uniform 
boundedness on B/\ x x B /\ 2 provided that it holds uniformly in Ai and A 2 ( i.e provided that 7 and /a 
in Definition 0 are independent of Ai and A 2 ). 

Remark 4 An interesting corollary of Theorem [^ is that the cascade of two USAS systems remains 
USAS, roughly, provided that the assumptions of Theorem Q hold with <5i = 62 = 0 and that 0 is 
replaced by the stronger condition: 

dVi dV\ 

— {t, x i) + — (t, x i)fi(1;, x i,di) < -k Al Vi(t, x i) , V |xi| > <5i, (11) 

A similar adaptation of Proposition 2 may be obtained. The proof of this statement is omitted for 
lack of space, but it follows along the same lines 3 as the proof of Theorem ^1 Furthermore, it is worth 
pointing out that, for an autonomous x \—subsystem (i.e., f\(x\,6\)), the original requirement 0 ) 
remains sufficient. 

We present a result that may help to check Assumption 0 in some particular contexts. It requires 
the non-positivity of the derivative of a Lyapunov function on a sufficiently large domain. The proof is 
given in Section mil 

Proposition 3 Let b be a positive constant. Suppose that there exists a continuously differentiable 
function V and two class /Coo functions a and ci such that, for all t £ R>o and all x £ R n , 

a(\x\) <V(t,x) <a(\x\) (12) 


x £ 7i(a, b) =k —(t,x) + —(t,x)f(t,x)< 0, (13) 

where a designates a positive number such that a(a ) < a(b). Then, for all to £ R>o, the solutions of 
0 satisfy 

l^ol A o' -1 ° a(b) =k \x(t,to,xo)\ < b, Vf > to . 


4 Application in robot control 


To illustrate the utility of our theorems, we consider the common problem of set-point control of 
a rigid-joint robot manipulator under PID control and and taking into account the dynamics of the 
actuators. The Lagrangian dynamics of a robot manipulator with n rigid-joints is given by 

D(q)q + C(q, q)q + g(q) = u (14) 

where D(q) £ W ixn is symmetric positive definite for all q £ R n , N(q,q ) := D(q ) — 2 C(q,q) is skew- 
symmetric for all (q, q) £ R n X R n , u £ R” corresponds to the input torques. As most common in the 
literature of robot control, we restrict our attention to systems satisfying the following (cf. 123 ED- 


Standing assumption 1 The functions D(-), C(-,-), g(-) are at least twice continuously differentiable 
and the partial derivatives of their elements are over-bounded by non decreasing functions of 4 |g| 
and \q\. Furthermore, we assume that there exist positive constants d m , dM and k c such that 5 for 
all q and q £ R n , 


d m < \D(q)\ < d M , \C(q,q)\ <k c \q\ , 


dg(q ) 
dq 


— kg ■ 


3 The application of Lemma Q (see Section If). 1 1 is no longer required in view of 0J- 

4 Notice that this is true for matrices containing polynomial and trigonometric functions which is fairly common in 
Lagrangian models of physical systems other than mechanical -cf. |3inr|. 

5 This is true for instance for open kinematic chains with only revolute or only prismatic joints. See e.g. E 21 IIH] 






We consider that the input torques u £ K" are delivered by Direct-Current (DC) motors, whose 
dynamics are given by 

Li + Ri + k b q = v (15) 

where i € M n is the vector of rotor currents, L and R are the rotors’ inductances and resistances 
respectively 6 , hqj with j < n is the back-emf voltage in each motor and v is the vector of input 
voltages, i.e., the control inputs. Each motor produces an output torque Uj = k t ij with kt > 0. 

Our control problem is to design v so that the robot manipulator stabilizes at a desired constant 
set-point (q = q*, q = 0). The aim is that the robot coordinates approach the reference operating point 
from any initial conditions in an arbitrarily large set (limited only by physical constraints). Furthermore, 
it is imposed that control be of the PID type. Accordingly, disregarding the DC motor dynamics, the 
input torques that achieve the control objective are given by 

u* = —k p q — k d q + v (16a) 

z> = ~hq, zz(0) := g{q*) (16b) 

where k p , k d and ki are positive design control gains and g(q*) is the best “guess” of the unknown 
constant pre-computed gravitational forces vector. 

We stress that the above setting is fairly common in practice of robot control; not only PID control is 
probably the most popular control technique but, often, industrial manipulators come with a black-box 
controller of PID type, meaning that control design for the user of an industrial robot boils down to 
gain-tuning for the built-in PID. 

Here, our control objective is achieved via cascaded-based control, relying on our results on USPAS. 
The approach consists in designing a reference i* := u*/kt (so that, when i := i — i* = 0, we have that 
u = u*) and building a control law v that makes i go to zero. 

Proposition 4 Under Standing assumption Hi the system (HU), hu in closed-loop with HU and 

u* 

v : = R'i + Ri* + k b q + Li* , i* = — . 

kt 

is uniformly semiglobally asymptotically stable (USAS). □ 

Proposition HI establishes that, if one knows how to semiglobally asymptotically stabilize a robot using 
PID control when neglecting the DC drive dynamics, then the same stability property can be established 
when these dynamics are taken into account. In other words, we claim that, given any domain of initial 
errors, one can always find control gains (namely k p , ki , k d ) such that the point (q = q*,q = 0) is 
uniformly asymptotically stable on this set of initial conditions. Moreover, the tuning of the PID 
control gains can be made disregarding the DC drive dynamics. 

Sketch of proof. (Proposition ^ For analytical purposes, let £\ > 0 be sufficiently small and define the 
variables s := q + jr(</((?*) — v) and k p \= k p — > 0. Then, the closed-loop system can be written 

in the following cascaded form: 

D(q)q + C{q , q)q + g{q) - g{q*) + k' p q + k d q - k t s 

s 

i 

Then, the proof of the proposition can be constructed by applying Theorem HI and using Remark HJ 
For this, three basic properties must be shown: 1) the motor closed-loop system (1181) is USAS, 2) the 

6 For simplicity and without loss of generality we consider that we have n identical motors -i.e. same resistance, 
inductance, torque constant, etc. 


hi 

(17a) 

1 

q-\ - q 

(17b) 

£l 

RPR' . . 

r *’ 

(18) 
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robot system in closed loop with the PID controller, i.e. G3> with i = 0, is USAS and 3) the trajectories 
of the cascaded system are uniformly bounded. 

Sketch of proof of 1): This follows directly by observing that (1181) is uniformly globally exponentially 
stable for any positive R'. 

Sketch of proof of 2): This follows after lengthy calculations that we do not include here for lack of 
space and because they follow along similar prooflines as those in ESI. For the purpose of illustrating 
the use of Theorem Q we claim that the result can be established using the Lyapunov function 

Vi ■= tU 7 D(q)q + ^k' p \q\ 2 + U(q) ~U{q*) - q T g(q*) + ^l s 2 + £l q T D(q)q + £ 2 s T D(q)q 

where U(-) corresponds to the gravitational potential energy function and £2 > 0. We emphasize that, 
except for the last term and a slightly different notation, the function above comes from m and, more 
indirectly, from energy-like Lyapunov functions widely used in the robot-control literature. In particular, 
conditions for positive definiteness 7 and radial unboundedness hold under the Standing assumption 1 
and can be derived following m and some of the references therein. 

It can also be shown that, for any A > 0, there exist k p , kd and ki such that the time derivative of 
V\ along the trajectories of m satisfies, for all |xi| := \/\q\ 2 + |<?| 2 + |s| 2 < Ai and all t > 0, 


V < -y I9I 2 



£2 ki 2 

~ir s 


provided that £\ and £2 are chosen small enough, and that the choice of the gains can be made in the 
following manner 8 : 

kd — &d bd^ 1 , kp — a p -{- 6 pAi, ki — ai T 6 ^A^, ( 19 ) 

where ad, a p , ai, bd, b p , bi are positive constants. In addition, the bound on the gradient © directly 
follows from the smoothness of V\ and its time-independency. This establishes Assumption [21 when 
considering di = 0. 

Sketch of proof of 3): This can be established following the same analysis as above, with the same 
Lyapunov function, but considering the interconnection term kfi on the right hand side of the closed 
loop equation (117 al) . One obtains that, in this case, 


- 2 |y| 2 


e 2 ki s 2 


+ £1 \q\ + £2 |s|) k t \i\ 


which, in view of the uniform boundedness of i(t), satisfies V < 0 for “large” values of |xi|. This 
observation, combined with the linear dependency of the gains m in Ai, allows to fulfill Assumption 
01 in view of Proposition 01 


Remark 5 Cascaded-based control of robots taking account of the robot dynamics was first used in DU 
in trajectory control of manipulators with AC drives. In that problem, the motor dynamics is highly 
nonlinear and global asymptotic stability of the closed-loop system of m with the corresponding ideal 
control input u* is obtained. However, for the problem that we address here, we are not aware of any 
proof of global asymptotic stability of the closed-loop with PID control hence global results for 
cascaded systems fail in this setting. See also [l] for a result on control of robots taking into account 
the DC motors’ dynamics, but with knowledge of the gravity terms g{q). □ 

7 Positive definiteness of V may be established if kp > k g and if £1, £2 are sufficiently small. 

8 We can actually show that the V-set (cf. Definition 0 J of 11 71 is 

P( 0 , Ai) = {(kd, k p , ki) G R 3 : kd > ad + 6dAi, k' p > a p + b p Ai, ki > ai + LAi} ■ 
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5 Proofs 


5.1 Proof of Theorems |7| 

We start by introducing the following result, which is a direct adaptation of El Proposition 13] and 
allows V\ to be transformed into a more convenient form. 

Lemma 1 Let <5 be a nonnegative constant and X be a subset of M n \ B§. Suppose that there exist a 
continuously differentiable function V : R>o x X —■» R>o and some class ACoo functions a, a, a such 
that, for all x € X and all t > 0, 

a(|z| 5 ) < V(t,x) < a{\x\) 
dV dV 

X) + — (t, x)f{t, x) < -a(|x|). 

Then, for any positive k, there exists a continuously differentiable function V : M>o x X — * M>o and 
class /Coo functions a, a such that, for all x € X and all t > 0, 

Ql{\x\ 5 ) < V(t, x ) < a(|x|) (20) 


dV dV 


and, for any s € M>o, it holds that 


( 21 ) 


a 1 o a(s) = a 1 o a(s). 


“I „ -77“ ( 


( 22 ) 


If, in addition, there exists a continuous nondecreasing function c : M>o —> R>o such that, for all x £ X 
and all t > 0, 

dV 




< c(|*|), 


then there exists a continuous nondecreasing function c : M>o —>• M>o such that, for all x € X and all 
t > 0, 

dV 




< c x). 


(23) 

□ 


Proof. Following the prooflines of El Proposition 13], we see that the function V can be defined as 
poV where 

/ P( s ) = ex P (fi 5§j) = Vs > 0 
1 p{ 0 ) = 0 , 

and a is a convenient class /C function. The bound m can be established following the same reasoning 
as in the proof of El Proposition 13]. Furthermore, (1201) can be satisfied with a := poa and a := po a 
as p € ACoo, and we therefore have that 

oT l o a(s) = (p o a) -1 o (p o a) (s ) = (a -1 o p -1 ) o(poa) ( s ) = a -1 o a(s). 


Concerning the bound on the gradient we have that, for all x € X and all t € M>o, 


, 


< 


2V(a) 

a(V(x)) 


dV. . 


2a(|x|) .| -/i i\ 

< 7 n | U c( g ) < c a , 


where c(s) := ^[f)) c ( s )i which establishes the result. ■ 

Consider the function V\ generated by Assumption^ and let Lemma[T]with x = Ti(A.\,5i) generate a 
function Vi, class ACoo functions a^ x , ai^Ai, a positive constant ks lt a x and a continuous nondecreasing 















function c$ lt a : such that, for all x G Ax) and all t G M> 0 , (t2T)|) - (|2M|) hold for Vi- From © and 

©, it cilso holds that, for any 0, 

oJihi ° ^SiAi (^i) = 0 • ( 24 ) 

In the sequel, in order to lighten the notations, we refer to A| as simply a 1; a^Ax as ah, Ai as 
fei and Qi .Ai as ci. Even though no longer explicit with these new notations, the dependency of these 
functions in di and Ai should be kept in mind. For any given positive 5i, Ai, <5 2 and A 2 satisfying 
Ai > max{5i, Ao} and A 2 > 82 , let 7 Ai,a 2 be generated by Assumption 0 and define 

A := min {Ai; A 2 ; 7 Aj,a 2 (Ai, A 2 )} . (25) 

Next, choose any 0* G @1 satisfying Assumption [21 and any 0£ G D/ 2 (d 2 , A 2 ) n 0 2 given by Assumption 
01 We show that, provided that A, <5 2 are sufficiently small and that Ai, A 2 are large enough, there 
exists 8 G (0; A) such that Bs is UAS on B/\ for the system © with 9* = (0i,0 2 ). To that end, we 
first show that there exists a positive 83 such that the ball Bs 3 is uniformly stable. More precisely, we 
construct 77 G /Coo and 83 > 0 such that, for all xo G £> a , 

\xi(t,t 0 ,x 0 ,6*)\ S3 <rj(\x 0 \). (26) 

Then, we use this property to prove that a ball, larger than Bs 3 , is UA on / 3 a and we construct a ICC 
estimate for the solutions. Finally, we show that the estimates of the domain of attraction and of the 
ball to which solutions converge can be arbitrarily enlarged and diminished respectively. 


5.1.1 Proof of uniform stability of a ball 

The total time derivative of V along the trajectories of 0 with 9 = 9* yields 

V, = + ^(/i(i,zi,0*) + g(t,x,9*)x 2 ). 

Therefore it holds that, for all x\ G 'H{ 8 \, Ai) and all t > 0 


l>i < 
< 


-fciVi + 


avi 


\g(t,x,9*)\ \x 2 \ 


dx] 

-hVi +ci(|xi|)G(|x|) |x 2 | 


Defining 

F := {t>t 0 \ 81 < \xi(t,t 0 ,x 0 ,9*)\ < Ai} , 

and using the shorthand notation x\(t) for xi(t,to,xo,9*) and v\(t) := Vi(i, xi(i)) we get that, for any 
xo G Ba and any f G T, 

vi(t) < -kivi(t) + ci(|xi(t)|)G(|x(t)|) |x 2 (t)| . 


From Assumption 0 and in view of (ESJ, we can see that, for all xq G /3a, it holds that x(t) G Z3ax • 
Hence, for all t G T, 

vi(t) < -k\vi(t) + ci(Ai)G(Ai) |x 2 (f)| . 

Let Assumption [T] generate a class ICC function /? 2 such that 9 for any x 2 o G / 3 a 2 and any t > to, 


X 2 (t)\ < P 2 (12-201 ,t - t 0 ) + 82 ■ 


It follows that, for all xq G £>a and all i G I, 

vi(t) < -kivi(t) + ci(Ai)G(Ai)[/3 2 (|x 20 | ,t-t 0 ) + <5 2 ] 

9 Notice that |s| a < b •<=> |s| < a + 6, Vs £ R", a, 6 > 0 . 


( 27 ) 








which implies that 


x{t) £ H(Si, Ai) => hi(t) < -kivi(t) + c 3 (|x 0 |), 


(28) 


with 

C 3 (s) := ci(Ai)G(Ai)(/? 2 (s, 0) + 6 2 ), Vs > 0 . 

The rest of the proof of uniform stability consists in integrating £U over r, in order to construct a 
bound like m- To that end, we introduce the following tool, which may be viewed as a comparison 
theorem for differential inequalities that hold only out of a ball centered at zero. 

Lemma 2 Let 8 be a nonnegative constant and A be a subset of W 1 containing B$. Assume that there 
exists a continuously differentiable function V : M>o x M n —> M>o, class /C 00 functions a and a, a positive 
constant k and nonnegative constant c such that, for all x £ X and all t £ M>o, 

a(M 5 ) < V(t,x) < a(\x\) 


and, for all xo £ M n and all to G R>o ; 

x(t,t 0 ,x 0 ) £ X\B S => V(t,x(t,t 0 ,x 0 )) < -kV(t,x(t,t 0 ,x 0 )) + c. 

Then, for all xq £ R n and to £ M>o such that x(t, to, xo) G X Vt > to, we have that 

\x(t,t 0 ,x 0 )\ g < aT 1 (a(8) + D +a ~ 1 ^a(|x 0 |)e _fe(t_ * o) + , Vt > t 0 . 

□ 

Proof. For simplicity, we write x(-,to,xo) as x(-) and we define v(-) := V(-,x(-)). We distinguish two 
cases: whether the trajectories start from outside or inside B 5 . 

Case 1 : |xo| > 8 . 

In this case, there exists 10 To £ (0;oo] such that |x(t)| > 8 for all [to;to + To) and |x(to + To)| = 8 . 
Hence, using the comparison lemma, we get that v(t) < (v(t 0 ) — | for all t £ [to; to + To). 

Using the bounds on V, it follows that 

|®(*)| 5 < aT 1 (a([x 0 |)e _fe(t_to) + , Vt £ [t 0 ; t 0 + T 0 ). 

In addition, for each t > to + To, either |x(t)[ < 8 in which case 11 |x(t)| 5 < or 1 (a(d) + c/k), or 
|x(t)| > 8 . In this second case, we can again invoke the continuity of the solution to see that there 
exists a nonempty time-interval [r; r + T], with T £ (0; 00 ], containing t and such that |x(s)| > 8 for all 
s £ (t; t T T], with |x(r)| = 8 . Hence, integrating from r to t £ [r; r + T], we obtain in the same way 
as before that, whenever |x(f)| > 8 , it holds that 

\x(t)\ s < aT 1 {a( 8 )e~ k ^ + £) < a" 1 (a(\x 0 \)e~ k ^ + . (29) 

To sum up, for all t > to, we have the following: 

|x 0 | >8 => |x(i)| 6 < a -1 (a{ |x 0 |)e _fc(i ~ to) + ^ . (30) 

Case 2 : |xo| < 5. 

In this case, as long as |x(t)| < <5, we triviallj 5 ^ have that \x(t)\ s < a -1 (a(d) + c/k). If |x(t)[ > <5 at 
some instant t > to, then, again, there exists a nonempty time-interval [r;r + T], with T £ (0; 00 ] and 
r > to, containing t and such that |x(s)| > 8 for all s £ (r; r + T], with |x(r)| = 8 . Thus, from (I29f) . we 
obtain that 

l*(t)l« < « _1 ( a( 8 )e~ k{t - T) + < 0 C 1 ( a( 8 ) + . 

10 If |a:(t)| > 5 forever after, we consider that To = °o- 

11 This is direct by noticing that a(s) < a(s) for all s € K>o and that c/k > 0. 





Hence, for all t > to, 

|x 0 |<(5 =>- |x(t)| 5 < of 1 (a(5) + . (31) 

The conclusion follows from PH) and m ■ ■ 

Applying Lemma [21 to (E%1) with V = V\, k = k\. c = c 3 (|xo|) and A” = Ha, we get in view of 
Assumption 0 and (E51) that, for all xo £ Ha and all to £ R>o, 

\x(t)\ Si <a f 1 + C3 ^°^ +«r X ^«i(kol) + C3 ^°^ ’ Vt>t 0 . 

Define the following: 

Jf _ £ , —1 ( Z' £ \ | C 3 ( 0 )^ , -1 ( £ 3 ( 0 )^ 

<5 3 <3.+a, \a 1 (6 1 ) + —)+a 1 J 


— <3i + 1 (<5i) + 


ci(A 1 )G(A 1 )<y 2 ^ _! fc^A^A^ 

' —X 


Ah 


Ah 


(32) 


and, for all s £ M>o, 


^ —1 (— /£ \ , c 3( s )\ , —1 (— f \ i c 300 > \ -1 (- ts ^ , c 3(0)^ _i Z'c 3 (0)\ 

r?(s) -“- 1 + (^i(s) + -^J -oh \a 1 (S 1 ) + -j-a 1 j 


We then conclude that, for any xo £ Ha and all to £ M>o, it holds that 

|®i(*)li 3 < »7(|®o|), Vt>t 0 . 


(33) 


Uniform stability of B § 3 on Ha follows by noticing that rj is a class K. function. This can be seen by 
recalling that c 3 is a continuous increasing function. 

5.1.2 Proof of uniform attractiveness of a ball 

Consider again P7|) . Since is a KC function there is a time t 3 > 0, independent of to and xo, 
such that 

/? 2 (A, t — to) < S 2 , Vt>to + ti- 
Hence (1771) implies that, for all f £ T fl R> to+tl and all xo € Ha, 

i»i(t) < -feixi(t) + 2 ci(Ai)G(Ai)5 2 • 

Applying again Lemma [21 and recalling that, from Assumption EJ |xi(to + ti)| < Ai, it follows that, for 
all xo £ Ha, all to £ M>o and all t > to + ti, 

|x(t)| 5l < ar 1 (a^i) + 2Cl(Al) fc ^ (Al)<52 ) + ^ ^Qi(|g(*o + ti)\)e~ kl ^~ to ~ tl] + 

< «r X + 2Cl(Al) fe G(Al ^ 2 ) + a" 1 ^ai(A 1 )e- fcl ( t - t °- t i) + 2ci(Al )G(A 1 )5 2 ^ 


Defining 


1 ( ai(Ai) 

h ■= ti + -r- In _ ( r r 
«i V cki(oi) 


we then see that, for all xo £ Ha, 

|xi (t)| < 64 := Si + 2 off 1 ^«i(< 5 i) + 

In other words, we have that 

|xi(f,f 0 ,xio)| 54 = 0, Vt>t 0 + t 2 


2c 1 (A 1 )G(A 1 )<y 2 

Ah 


Vt > to +1 2 • 


(34) 
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Finally, let 


6 := max{J2 ; £3 ; £4} 


(35) 


Then we see that (13311 implies that |xi(t )| 5 < r;(|xo|) for all t > to- From this and what precedes it is 
not hard to see that, for all xq E £>a, 


xi(t)\ s < v(\xo\)e {t to tz) , Vt>t 0 - 


Thus, recalling that t 2 depends neither on to nor on xq, and defining 


/3(s, t) := max |ry(s)e ^ ; ^(s, t)| , Vs,f> 0 


= max 


we conclude that, for all xq € 13a, 


x{t)\ s < P(\x 0 \,t - t 0 ), Vt>t 0 . 


UAS of £>5 on 13 a follows by noticing that /? is a class /C£ function. 


5.1.5 Semiglobality and practicality 

It is only left to show that 6 and A can be arbitrarily reduced and enlarged respectively. It follows 


directly from m and 031) that, by picking Ai and A 2 large enough, A can be made arbitrarily large. 
Concerning 5, we see with (1331) and (1331) that 



Hence, in view of (1331) and recalling that ci, G, k\ and a t are independent of 62 , we see that, for the 
chosen Ai and A 2 , both <5 3 and 64 can be taken as small as wanted by picking 5] and 62 sufficiently 
small. Hence, in view of (1331) . it is also the case for 6 . 

Thus, it suffices to pick the parameters 0* and 0\ generated by the chosen <5i, Ai, 62 and A 2 , to 
conclude that, for any A > 5 > 0, there exists some parameters 0\ E ©i and 0\ £ ©2 such that B$ is 
UAS on Ha for the system m with 0 = 0 *, which establishes the result. 

5.2 Proof of Proposition 0 

We prove this result by showing that it actually constitutes a special case of Theorem [l] To this 
end, first notice that applying Lemma Q (see Section [5. II) to V ensures the existence of a continuously 
differentiable function V such that, for all x E 7V(<5, A) and all t E M>o, 


«*,a(I®I*) < V(t,x) < ai,A(|x|) 


< -4,AV(t,x) 

dV, , _ „ 1N 


(36) 


— (t,x) < c 5i a(|x|) 


hold with a positive kg t a, a continuous nondecreasing function c^a and some class K , 00 functions A 
and as,A satisfying 

afX 0 «< 5 , a ( s ) = o as,A{s) , Vs > 0. 

Inverting the two sides of this inequality yields: 


a s ,A 0 A (s) = 0 a (s), Vs > 0 . 


13 





Consequently, in view of © and ©, we have that, for all A > 0, 

liiri aji oaj AM) = 0 
<5^0 —l 


i-l 


and, for all <5 > 0, 

Based on this, let A be any given positive constant and choose 5 small enough that 


lim a SA oas A (A) = oo . 

A—xx) 


( 37 ) 


“ 5 ,A 0 < A - 

Then, the requirements of Proposition [3] are fulfilled and we get that 

|*o| < ° « 5 ,a( a ) => \ x (t)\ A A , Vt>t 0 . 

In view of (EH) , we then see that the solutions of © satisfy a uniform boundedness as in Assumption 
Moreover, from m , and defining v(t) := V(t, x(t, to,xo, 9)) and T := {i > #o I x (t, to,xo, 6) G Ti.(5, A)}, 
we get that 

v(t) < —ks tA v(t) , Vi G T, 

which corresponds to ®) with C 3 (s) = 0. The rest of the proof follows along the same lines as in 
Section Oby noticing that both Assumption [21 and the bound on the gradient were only used in order 
to establish C9> . 


5.3 Proof of Proposition [3| 

We claim that, whenever V(t,x) = a(b ), its derivative along the trajectories of ©, which we 
denote by V, is non positive. To this end, notice that m implies that, if V(t,x) = a(b), then 
x G 7i (a -1 o a(b), b), which is nonempty (since a(b) < a(b)) and included in Tt(a, b ) (since it is assumed 
that a (a) < a(b)). Hence, the claim is proved in view of (THT1) . For any 1 0 G M>o and any xq G M”, by 
defining v(t) := V(t,x(t,to,xo)), we therefore get that, for all t > to, 

v(t) = a(b ) => v(t) < 0 , 

which ensures in its turn, by the continuity of v(t), that 

v(t 0 ) < a(b) =>■ v(t) < a(b), Vt>to- 
The conclusion follows by noticing that, from (El, 

|*o| < oT 1 0 Qfb) =>• v(to) < a(b) , and v(t) < a(b) |x(t)|<6. 

6 Conclusion 

We have presented results for uniform semiglobal practical asymptotic stability of nonlinear time- 
varying systems. Our main theorem establishes that the cascade of two USPAS systems remains USPAS; 
it relies on a condition of local boundedness of the solutions of the cascade and the knowledge of a Lya¬ 
punov function for the perturbed subsystem, in the absence of the interconnection term. As a corollary, 
we have that, under similar conditions, the cascade of two uniformly semiglobally asymptotically stable 
systems remains USAS. The latter generalizes, in its turn, previous results reported in the literature. 
Further research is carried out in the direction of extending these results to global practical properties 
and on deriving sufficient conditions for boundedness of solutions. 

We have also illustrated the usefulness of our main results in the PID control of manipulators with 
external disturbances. The case-study presented here addresses, as far as we know, an important open 
problem in the literature of robot control. 
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